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ABSTRACT 

Faraday rotation measures (RMs) of the polarized emission from extragalactic radio sources oc- 
culted by the coronal plasma were used to infer the radial profile of the inner heliospheric mag- 
netic field near solar minimum activity. By inverting LASCO/SOHO polarized brightness {pB) 
data taken during the days of observations on May 1997, we retrieved the electron density dis- 
tribution along the lines of sight to the sources, thus allowing to disentangle the two plasma 
properties that contribute to the observed RMs. By comparing the observed RM values to those 
theoretically predicted by a power-law model of the radial component of the coronal magnetic 
field, using a best-fitting procedure, we found that the radial component of the inner heliospheric 
magnetic field can be nicely approximated by a power-law of the form S, = 3.76 G in a 

range of heights from about 5 to 14 Rq. Finally, our analysis suggests that the radial computation 
of the potential field source surface (PFSS) model from the Wilcox Solar Observatory (WSO), 
assuming a radial field in the photosphere and a source surface located at Rss = 2.5 Rq, is the 
preferred choice near solar minimum. 
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1. Introduction 

Despite its fundamental importance, no direct information is available as yet on the inner helio- 
spheric magnetic field, the innermost in situ measurements still being those obtained in the seven- 
ties by the two Helios probes at ~ 62 Rq. Such crucial measurements, in the region where the solar 
wind is heated, accelerated and decoupled from the coronal plasma, will be only available with the 
advent of the Solar Probe Plus missiorj^ planned by NASA to be launched no later than 2018. On 
the final three orbits. Solar Probe Plus will be able to fly within 8.5 Rq from the Sun's surface. 
Empirical estimates of magnetic fields in the inner heliosphere are only possible very close to the 
base of the corona (e.g., Lin et al. 2000; Lee 2007), but their outward extrapolation often involve 
simpUstic potential (or force-free) assumptions and the hypothesis of low plasma which may 
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not be appropriate in the outer corona (Gary 2001). Under appropriate assumptions, magnetic field 
strength estimates in the outer corona have been obtained from the analysis of the radiation emitted 
during the passage of coronal shock waves (e.g., Mancuso et al. 2003; Vrsnak et al. 2004; Cho et 
al. 2007; Bemporad & Mancuso 2010, Gopalswamy & Yashiro 2011; Kim et al. 2012). However, 
one of the best observational methods for obtaining information on the strength and radial profile 
of the magnetic field in the inner heliosphere still remains the analysis of Faraday rotation measure- 
ments of linearly polarized radio signals from galactic or extragalactic radio sources (e.g., Sakurai 
& Spangler 1994; Mancuso & Spangler 1999, 2000; Spangler 2005; Ingleby et al. 2007; Mancuso 
& GarzelU 2007; Ord et al. 2007) or from the transmitter of a spacecraft (e.g., Stelzried et al. 1970; 
Patzold et al. 1987; Jensen et al. 2005). For a fairly recent comprehensive review of coronal Faraday 
rotation observations, see Bird (2007). 

In order to understand the physics behind the Faraday effect, it is useful to consider a linearly 
polarized wave as two counter-rotating circularly-polarized waves. Faraday rotation occurs because 
the phase velocity of the left circularly polarized component travels faster along the magnetic field 
than the right circularly polarized component, thus resulting in a net rotation of the wave's 
polarization position angle given (in radians) by: 



In the above equation, expressed in cgs units, A is the wavelength of the radio signal, is the 
electron density, B is the vector magnetic field, ds is the vector incremental path defined to be pos- 
itive towards the observer, e and rrie are the charge and mass of the electron, and c is the speed of 
light. The above expression can also be written as = ??RM, where RM is the rotation measure 
(RM), defined to be positive for B pointed towards the observer. Essentially, RM yields informa- 
tion on the integrated product of the hne of sight (LOS) component of the magnetic field and the 
electron density. Since Faraday rotation is an integrated quantity and proportional to the product 
of two independent quantities, in order for the magnetic field component of the observed RM to 
be determined, there must be an independent determination of the electron density, whether from 
observations or models. For a spherically symmetric corona, that is, with a perfectly symmetric 
electron density and magnetic field distributions, the obvious result would be RM at all lati- 
tudes. This is due to the fact that the LOS component of the magnetic field would reverse at the 
point of closest approach to the Sun. Overall, the analysis of the Faraday rotation observations 
reduced by Mancuso & Spangler (2000) showed values different from zero, reaching a maximum 
value (in absolute magnitude) of -10.6 rad mr^ for one of the radio sources at ~ 6 Rq. According 
to the study of Mancuso & Spangler (2000) , this effect was mostly attributed to a non-symmetric 
distribution of the electron density along the LOS, thus working as a "weighting function" for an 
otherwise symmetric magnetic field distribution. Although the electron density distribution used in 
Mancuso & Spangler (2000) was not synnmetric, it was still an analytical approximation inspired 
by the work of Guhathakurta et al. (1996), with two different analytical expressions for the streamer 
and coronal hole components. In fact, because of this assumption, no attempt was made to actually 
derive, by a best-fitting procedure, an analytical expression for the radial component of the mag- 
netic field. Moreover, in that work, the position of the magnetic neutral line, which coincides with 
the heliospheric current sheet, was deduced from a potential field expansion of the photospheric 
magnetic data from the Wilcox Solar Observatory (WSO) and using a potential field source surface 
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(PFSS) model (Altschuler & Newkirk 1969; Schatten et al. 1969) with a source surface at a fixed 
height of 2.5 Rq. Although Mancuso & Spangler (2000) could obtain a fair agreement between 
RM model and observations, their analysis was only able to support the validity of the Patzold et 
al. (1987) coronal magnetic field model, which was used in input as an approximate expression for 
the radial magnetic field profile, at the heights covered by the RM observations. 

The aim of this work is to improve the previous analysis of Mancuso & Spangler (2000) by 
introducing direct information on the electron density distribution along the LOS for the above 
observations, thus allowing to obtain an estimate of the strength and radial profile of the in- 
ner heUospheric magnetic field. In this work, instead of using an analytical model for the elec- 
tron density distribution along the LOS, we will obtain this quantity empirically through inver- 
sion of polarized brightness (pB) measurements obtained from observations of the Large Angle 
and Spectrometric COronagraph (LASCO; Brueckner et al. 1995) telescope aboard the Solar and 
HeUospheric Observatory (SOHO; Domingo et al. 1995) spacecraft during the same days of ob- 
servations of the radio sources. In this way, we will be able to infer, by a best-fitting procedure, 
the radial profile of the magnetic field of the inner heliosphere that best accounts for the observed 
RMs. As a by-product of our study, we will also be able to evaluate the most suitable PFSS model 
among the three available models from WSO. 

This paper is organized as follows. In Sect. 2, we present some details of the observations and 
of the data reduction procedure. In Sect. 3, we introduce the magnetic field model. In Sect. 4, we 
compare the model with observations and discuss our results. Finally a summary and conclusions 
are given in Sect. 5. 

2. Observations and data reduction 

2.1. Radio observations 

The observations of the radio sources discussed in this paper were made on four days in 1997 
(May 6, 11, 22 and 26) with the Very Large Array (VLA) radio telescope of the National Radio 
Astronomy Observatory (NRAO). Each source was reobserved far from the Sun at the same fre- 
quencies to obtain the intrinsic RM and that due to the interstellar medium but without a coronal 
contribution. The frequencies used were 1465 MHz with a bandwidth of 50 MHz and 1665 MHz 
with a bandwidth of 25 MHz and were chosen to have sufficient separation in to allow for an ac- 
curate determination of the RMs. Ionospheric calibration was carried out to correct for ionospheric 
Faraday rotation. Two major factors limiting the accuracy of VLA polarization measurements are 
the removal of ionospheric Faraday rotation and the correction for instrumental polarization. The 
task FARAD of the NRAO Astronomical Image Processing System (AIPS) was used for iono- 
spheric RM estimation and correction. This task calculates the appropriate Faraday rotation correc- 
tion using a phenomenological model of the ionospheric electron density. This model (Chiu 1975), 
takes as input the monthly average Zurich sunspot number to estimate the ionospheric electron 
column density. The adequacy of the ionospheric correction was tested by observing four sources 
far from the Sun, where coronal Faraday rotation is negligible. For these sources the standard de- 
viation of the rotation measure fluctuations around the mean was about 0.3 rad nT^ (see Mancuso 
& Spangler (2000) for details). This quantity can be considered as a rough estimate of the mea- 
surement errors introduced by residual ionospheric Faraday rotation or instrumental polarization 
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Fig. 1. The above four panels show the location of the radio sources relative to the Sun on each of 
the four days of observation. The bull's eye symbol indicates the position of the Sun and the dotted 
line is the ecliptic. The size of the plotted symbol is a rough indicator of the absolute magnitude of 
the RM. Red squares correspond to positive RMs and blue circles represent negative RMs. 



errors. In both coronal and reference observations, calibration of the instrumental polarization (and 
subsequent correction) was achieved by observing a calibrator source over a large range in par- 
allactic angle. The procedure used is described in Mancuso & Spangler (1999). The observations 
of the radio sources discussed in this paper sampled solar elongations from ~ 5 to 14 Rq at vari- 
ous heliographic latitudes through 13 different lines of sight (see Fig. 1), with values for the RMs 
ranging from -10.6 to -1-3.3 rad m"^. More specific information on the VLA observations used for 
this work and on the data reduction can be found in Mancuso & Spangler (2000) and will not be 
repeated here. 

2.2. White-light observations 

Observations of the extended corona are primarily obtained by white-light coronagraphs that are 
able to observe coronal structures which are highlighted by the photospheric radiation Thomson- 
scattered by electrons in the ionized coronal plasma. The observable polarized brightness (pB) is 
defined as the difference in intensity of radiation polarized tangential to the limb of the Sun and 
radiation polarized normal to the limb. This quantity is directly related to the coronal electron 
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density by 

pB{x,y)= r Tir)neir,e,(f>)ds, (2) 
Jlos 

where is the electron density, T(r) is the Thomson scattering fimction, and the integral is carried 
out along the LOS. In general, it is difficult to infer the density distribution in the corona and inner 
heliosphere from remote sensing techniques since the signals usually derive from the contribution 
of a superposition of different structures integrated along the LOS in the optically thin intervening 
plasma. Coronal mass ejections can also drastically affect the pB data as well as disturb the overall 
structure of the corona in a temporal scale of a few hours. Notwithstanding the above observational 
caveats, transient phenomena are much less common during solar minimum conditions. Moreover, 
the overall structure of the corona is observed to be axially symmetric by a high degree during 
the periods with lower solar activity due to the dominance of the dipolar component of the global 
magnetic field of the Sun. The electron density distribution can be thus estimated with a good 
degree of confidence from the observed pB using a technique developed by Van de Hulst (1950) 
that is particularly appropriate in the near-symmetric density distribution observed during solar 
minima conditions (e.g., Gibson et al. 1999). In this work, the inversion of the pB data, recorded 
daily around the Sun by the LASCO C2 telescope that observes the solar corona between about 2 
and 6 Rq, was done by fitting the observed radial profiles at steps of 3° of heliographic latitude. 
These pB radial profiles were then inverted with the above mentioned technique to yield electron 
density radial profiles below about 6 Rq. Finally, radial power fits of the form Heir) = Ar~" + Br~^ 
were performed in order to extend the profiles to 1 AU. Once the electron density maps were 
produced, the data were placed in a rectangular grid with observation times converted to Carrington 
longitudes and the Carrington maps were finally resampled by using a moving smooth filter and 
then a moving average, both applied to all pixels. Apart from the assumption of the corona being 
essentially stationary during the period corresponding to our observations, the above reconstruction 
of the global coronal electron density implicitly assumes that the observed Thomson scattered 
radiation is dominated by a relatively thin layer of scattering electrons centered on the point of 
closest approach to the Sun. This is clearly a zero-order approximation, since the actual pB should 
be considered as a weighted sum of contributions to the Thomson scattered radiation from electrons 
along all the line-of-sight. However, it is a fairly safe assimiption given that the electron density 
distribution in the corona falls off with distance as a high-exponent power-law. 

We finally mention that although there was also availabiUty of LASCO C3 measurements, 
probing the corona up to about 30 Rq, during the time interval under examination, these data were 
not used. In fact, apart from being noisier than the LASCO C2 data, it is known that the contribution 
from the F (dust) corona cannot be neglected above ~5 Rq (e.g., Koutchmy & Lamy 1985; Mann 
1992; Morrill et al. 2006) and the separation between the K and F coronae is not straightforward, 
thus making the pB-based inversion method more problematic and the determination of the electron 
density more difficult. Finally, in the range covered by the LASCO C3 field of view (about 4 to 32 
Re), it is statistically more probable, with respect to the observations obtained within the LASCO 
C2 field of view, to find coronal mass ejection (CME) signatures which would significantly distort 
the observed pB profiles at all heliolatitudes. 
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3. Magnetic field model 

The coronal magnetic field during solar minimum activity is dominated by the low-order magnetic 
multipoles, with the largest contribution coming from the monopole and dipole components. In this 
work, as in Patzold et al. (1987), the radial component of the global magnetic field Byadir) is thus 
assumed to have the functional form given by: 

= ^ + ^. (3) 

The dual power-law form of eq. (3) for is the linear combination of a dipolar field (oc r^^'), 
representing the dominant component of the global solar magnetic field at solar minimum and a 
monopolar, solar wind component (oc r"^) prevaihng at large distances from the Sun. Closer to the 
Sun, higher-order multipole components are present, but their contribution can be ignored at the 
heights relevant to these RM observations, since such fields fall off rapidly with height. Analysis 
of HeUos data (Mariani et al. 1979) showed that the absolute value of daily averages for the radial 
magnetic field component scaled as By oc between 0.3 and 1 AU near solar minimum, thus 
justifying the presence of the monopolar component in eq. (3). Moreover, the radial component of 
the heUospheric magnetic field, as detected by Ulysses during the same solar minimum as the one 
investigated in this work, was found to be remarkably constant with latitude, with an average value 
of \Br\ ~ 3.1 X 10-5 Q (Smith & Balogh, 1995; Balogh et al. 1995). In this paper, the value of 
representing the scale factor for the solar wind component (oc r"^) at great distances from the Sun, 
was thus fixed by this observational constraint, that is, B02 ~ 1.43 G X Rq^ for Br expressed in 
Gauss units and r in solar radii (Rq). 

HeUos observations also revealed that during solar minimum the slow solar wind of the streamer 
belt is restricted to a region of about ±20° around the heUospheric current sheet. A similar result 
has been obtained by the Ulysses probe (e.g., Wock et al. 1997) that detected a sharp transition 
in latitude from slow to fast solar wind. Thus, in the following, we will assimie for the sake of 
simpUcity that the slow wind (streamer belt) region is confined in a region within 20° above and 
below the heUospheric neutral line, with the fast wind occurring beyond this latitudinaUy bounded 
region. Another important parameter is the Alfven radius r^i, that is, the distance at which the solar 
wind becomes super- Alfvenic in the outer corona. Below i'a, the field is strong enough to control 
the plasma flow and cause the solar wind to corotate with the solar wind that has no azimuthal 
component in a corotating frame. Above r^, the plasma is released from corotation and the field 
wraps up following, approximately, an Archimedean spiral (Parker 1958). For r < ka, in a spherical 
coordinate system {r,d,(p) the azimuthal component of the magnetic field, B^, is null. Beyond r^i, 
B^ is different from zero, slowly growing with r as oc QQ(r- r^) cos 6/Usw, thus being a function 
of the coronal rotation rate Qq, the solar wind speed m^, and of heUolatitude 6. The value of 
remains uncertain as yet, but it has been estimated as ranging from ~ 10 Rq over the poles to ~ 30 
Rq over the equator (Scherer et al. 2001). In the following, will be modeled by a step function 
assuming a value of = 30 Rq in the slow wind region around the equator and = 10 Rq in the 
fast wind region towards the poles. 

Observations of the solar coronal rotation rate in solar cycle 23 have shown that the extended 
corona is more rigid than the photosphere (e.g., Giordano & Mancuso 2008; Mancuso & Giordano 
2011, 2012) but still sensibly latitude-dependent. In this work, the coronal rotation rate has been 
approximated by the functional form given by Chandra et al. (2010) for the soft X-ray corona 
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in 1997 from images obtained with the soft X-ray telescope (SXT) on board the Yohkoh solar 
observatory. As for the solar wind speed, the value for the slow component was chosen on the basis 
of the average speed of the slow wind at 1 AU as deduced from the data obtained by the Solar Wind 
Experiment (SWE; Ogilvie et al. 1995) on the Wind spacecraft during May 1997. As for the fast 
solar wind at 1 AU, we relied on the analysis of Quemerais et al (2007), who deduced the velocity 
profiles for Spring 1997 from the combination of SWAN, LASCO, and IPS data. We remark that 
although we took into account, for completeness, the contribution of the azimuthal component of 
the magnetic field to the observed RM, this amount remains practically negligible (of the order of 
a few percent) in the range of heights considered in this work. 

Finally, a necessary ingredient in coronal RM studies is the knowledge of the location of the 
heliospheric current sheet with respect to the LOS to the occulted radio sources. This informa- 
tion is important to infer the polarity of the magnetic field along the integration path to the different 
sources, especially near the solar equator. This is customarily done through the use of PFSS models. 
Because of their simplicity, PFSS models are still used extensively within the solar and heliospheric 
commimities. Riley et al. (2006) demonstrated that PFSS models are especially applicable during 
solar minimum activity and that can generate solutions that closely match those generated by MHD 
models for cases when time-dependent phenomena are negligible. PFSS models assume that the 
magnetic field can be described as the negative gradient of the scalar potential <I> which satisfies 
Laplace's equation V^O = 0, thus implying that V x B = Anijc - 0, so that current densities 
J are neglected. Laplace's equation for O is solved in the space between the photosphere and an 
outer spherical shell (the so-called source surface, with radius R^s), yielding a solution in terms 
of Legendre polynomials. The so-obtained scalar potential is uniquely determined given irmer and 
outer boimdary conditions at the photosphere and at the source surface, respectively. The magnetic 
neutral Une, when propagated outwards, can be used as a proxy for the latitudinal location of the 
warped heUospheric current sheet, and this location can be readily and reliably obtained through 
the PFSS models. The coronal magnetic field calculated from photospheric field observations with 
the PFSS models is tabulated by the WSO and available online at their website with two differ- 
ent extrapolation methods. The classic computation, which locates the source surface at 2.5 Rq, 
assumes that the photospheric field has a meridional component and requires a somewhat ad hoc 
polar field correction in order to better match the in situ observations at 1 AU. The radial computa- 
tion, with located at 2.5 Rq and 3.25 Rq, assumes that the field in the photosphere is radial and 
requires no polar field correction. The position of the magnetic neutral line, separating regions of 
opposite polarity, is simply given by the contour of zero magnetic field strength. In our model, the 
magnetic field is thus assumed to be almost purely radial above the source surface and a function 
only of radius except for the polarity switching sign across the neutral line. The location of the he- 
liospheric neutral sheet and, consequently, the polarity of the magnetic field, will be inferred from 
the position of the magnetic neutral line as deduced from the three available PFSS models from 
WSO. 

4. Results and discussion 

In the following, we will implement three sets of models. In the first set of models (Case A), in 
order to reproduce the observed RMs, we will use the functional form for the radial component of 



7 



S. Mancuso and M. V. Garzelli: Radial profile of the inner heUospheric magnetic field 



the heliospheric magnetic field given by eq. (3) with Bqi left as a free parameter. We will further 
assume that the radial variation of the electron density with heliographic coordinates can be 
reproduced with accuracy by the pB inversion method outlined in Sect 2.2. This hypothesis might 
be questionable in that it would require both perfect caUbration of the available LASCO C2 data 
and that the extrapolation to higher distances is accurate enough. In a second set of models (Case 
B), we will set Bqi = and investigate the possibility that the electron density n^. has been obtained 
only up to an unknown multiplicative factor C„, so that its magnitude needs some adjustment 
(that is, tie — > Cnne) to be in agreement with the actual RM observations (see also the discussion 
in Ingleby et al. 2007). By setting the parameter Bqi to zero and letting C„ to be the only free 
parameter, we will test the possibihty that the radial magnetic field, down to a few solar radii from 
the Sun, might be simply described by a solar wind component (oc r~^) (e.g., Vrsnak et al. 2004; 
Spangler 2005; Ingleby et al. 2007). Finally, in the most general case (Case C), both C„ and Bqi 
will be left as free parameters. All these models will be separately tested by assuming the three 
available PFSS computations from WSO (i.e., the classic computation with R^^ = 2.5 Rq and the 
radial computation with R^^ - 2.5 Rq or 3.25 Rq). 

The unknown parameters in the above described models will be evaluated by minimizing the 
weighted sum-of-squared-residuals: 



where RMobsj is the observed RM of the ;-th LOS, RMi is the model RM for the same LOS, and w, 
is a weight assigned to each measurement depending on its uncertainty. Actually, the proper choice 
of (Oi is a critical issue in our statistical analysis. In fact, the available uncertainties ctrmj - o'i,obs 
(of order ~ 0.1 - 0.3 rad m"^ for most observations) for RMobs,i as quoted in Mancuso & Spangler 
(2000) , were assumed to be solely due to radiometer noise. The RM fluctuations attributable to 
coronal turbulence and MHD waves, 6RM, can be, however, much higher than <trm,c For example, 
RM measurements from the Helios data showed RM fluctuations up to 6RM ~ 2 rad m"^ with a 
time scale of roughly half an hour (HoUweg et al. 1982; Bird et al. 1985). Similarly, Sakurai & 
Spangler (1994) set an upper limit of 5RM ~ 1.6 rad mT^ for a radio source observed during solar 
maximum whose LOS passed within ~ 9 of the Sun, although Mancuso & Spangler (1999), in 
an analogous analysis near solar minimum, found SRM < 0.4 rad mT^. The overall level of turbulent 
RM fluctuations is thus not well constrained as yet and might be dependent as well on the particular 
phase of the solar cycle. Since there was no simultaneous independent estimate of the contribution 
of RM fluctuations attributable to coronal turbulence and MHD waves during the days of obser- 
vations, the uncertainties £rj?M,/as quoted in Mancuso & Spangler (2000) , may underestimate by 
some factor the true uncertainties in each measurement, given by (Trm/ = (o"smi SRM^Y^^. As 
a consequence, the true uncertainties ctrm/ remain actually unknown. According to the previous 
discussion, in the process of minimization for the evaluation of the unknown parameters, we 
have considered different o), estimates, by proceeding through three difl'erent steps, as explained 
below. 

First of all, in Step 1, we analyzed the RM data by considering as their total uncertainties 
crmj the ones that can be ascribed to radiometric noise only, i.e., we completely neglected the 
possible contribution of coronal turbulence and MHD waves to RM fluctuations. We minimized 
the for each of the cases A, B and C, and for each available PFSS model. We then computed 




(4) 
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the reduced chi- squared = x^l^ by taking into account the number of degrees of freedom (v) 
in the different cases (12 for Cases A and B, and 11 for Case C, due to the different number 
of estimated parameters). We also estimated the goodness of fit by computing the probability of 
observing a minimum of the distribution that arises from eq. (4) by varying the model parameters, 
larger than the one actually observed, assuming that this minimum is actually distributed according 
to ?i with V degrees of freedom (e.g., Garzelli & Giunti 2002). Calling a this probability, that 
in practice represents the area of the right tail of the distribution in the interval between [x^i„ 
and +oo), we say that the fit is acceptable at 100% a confidence level. In all cases, the reduced 
chi-squared for the best-fit points were xl > 30, thus yielding a very poor a. These results clearly 
suggest that the uncertainty component due to radiometric noise largely underestimates the total 
uncertainty, at least for some of the sources, and that the uncertainties 6RM due to RM fluctuations 
cannot be neglected in our study. Secondly, in Step 2, in the absence of an experimental estimate of 
6RM for the measurements at hand, we considered the range of possible 6RM estimates provided 
in hterature (see the discussion above), in order to obtain a gross estimate of the typical order 
of magnitude of 6RM. Taking into account the interval spanned by these estimates, we decided to 
adopt an average value 6RM = 1 rad mT^, that we assumed to remain constant during the days of our 
RM observations. We also assumed that, as expected in a realistic situation, the RM fluctuations 
do not act by systematically increasing or by systematically decreasing all RM observed values 
in the same direction, but may act in a different way on each RM observation, i.e., they do not 
introduce any correlation between the RM observations of different sources. We thus replaced the 
uncertainties ct^m,/ = cTi^obs used in Step 1 with the uncertainties ctrm/ = ^^^obs SRM^y^^, 
flnding;^fy values compatible with 1, thus leading to a more meaningful statistical analysis. In this 
conflguration, we were able to identify the best PFSS model among those proposed by WSO (that 
turned out to be the same for all three cases A, B and C), the best-flt parameters for each case, and 
the uncertainties on them at 90%, 95% and 99% confldence levels. These uncertainties are given 
in the foUowing as confldence intervals in the cases where we just vary a single parameter (cases 
A and B), and as confldence regions when, instead, we vary more parameters at the same time 
(Case C). Finallly, in Step 3, we considered a statistical analysis with uniform weighting, i.e., with 
(jL>i arbitrarily set to 1, so to yield an unweighted least-squares. The results of this analysis turned 
out to be similar to those of the analysis of Step 2, thus proving the marginal role of cTi obs on the 
analysis' outcome, with respect to the dominant contribution due to turbulence and MHD waves 
(i.e, o-rm/ ~ SRM » ai^obs). 



In our opinion, the results obtained using the weighting criteria of Steps 2 and 3 are actually, 
from a physical point of view, more representative of the real situation, insofar the true imcertain- 
ties are most probably dominated by the contribution of RM fluctuations attributable to coronal 
turbulence and MHD waves. In the following, we will therefore concentrate our discussion mostly 
on the results concerning Steps 2 and 3. We remark, however, that the above conclusion is strictly 
valid only in the assumption (yet to be verified in the range of heliocentric distances considered in 
this work) of 6RM~ constant both in latitude and distance from the Sun. 
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Fig. 2. Scatter plot of predicted RMs versus observed RMs in the best-fit parameter configuration, 
using a model with Bqi left as free parameter and setting C„ = 1. Results are displayed for both 
cases of non-uniform weights, after the two analyses described in Step 1 (yellow squares) and Step 
2 (green circles), respectively. The dashed 1:1 line shows the diagonal, corresponding to the ideal 
situation where predicted and measured values are equal. The solid (red and brown) error bars de- 
note the o-RM.i uncertainties for each measurement due to radiometer noise as quoted in Mancuso 
& Spangler (2000), whereas the dashed (blue) errorbars denoted the ctrm/ total uncertainties ob- 
tained after adding the contribution of RM fluctuations, estimated to be 6RM - 1 rad m"^ for the 
sake of this analysis. 



4.1. Model with Bqi ^OandC„^l (Case A) 

To find the regions of parameter space consistent with the RM observations, we performed a first 
grid search allowing the two unknown parameters Bqi and Rss to vary within a plausible range 
and fixing the multiplicative factor C„ to unity (that is, assuming that the radial variation of the 
electron density with heliographic coordinates is reproduced with accuracy). The ranges of the 
model parameters that define our search space were constrained to be within a plausible interval of 
< Bqi < 20 (with Bqi expressed hereafter in units of G x Rq^) at Steps of 0.01 Rq. The result 
of the best-fitting procedure yielded a very well defined minimum, thus allowing an accurate 
determination of the best-fit parameters. The best-fit parameters to the observed data for Case A 
are Bqi - 10.54 and Bqi - 10.95 for the two weighing criteria explained above (Steps 2 and 3) and 
for the radial computation with Rss = 2.5 Rq. The radial computation with Rss - 3.25 Rq yields 
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Fig. 3. Same as Fig. 2, but using a model with C„ left as free parameters and setting Bqi - 0. 



higher values, while the classic computation with R^s = 2.5 Rq yields a very high value, 
hinting to Bqi - 0. How well a model performs compared to the measurements is conveniently 
assessed from the scatter plot of observations versus predictions. Tighter scatter of the data points 
around the 1:1 line obviously indicates more accurate model fit. In Fig. 2, we show a scatter plot 
of predicted RMs versus observed RMs for Case A obtained by using the best-fitting parameters 
and the two non-uniform weighting criteria listed in Steps 1 and 2. The estimated RM values 
look to be nicely distributed around the line 1:1 of perfect fit for all the range of measured RM 
values. The largest discrepancies between the results of the models of Steps 1 and 2, that differ 
because of the RM uncertainties attributed to each observation, also shown in the plot, are visible 
in the RM predictions for those sources that have a larger absolute experimental RM value. In 
particular, adopting the model of Step 2 appears to greatly improve the agreement between theory 
and experiment for the source with an highly negative measured RM (< -10 rad m"^). 



4.2. Model with Bqi = and C„ + 1 {Case B) 

In the model presented in the previous section (Case A), we assumed that the radial variation of 
the electron density with heliographic coordinates could be reproduced with accuracy by the pB 
inversion method applied to the white-light data, thus requiring perfect calibration of the available 
data and accuracy in the extrapolation of the LASCO C2 data to higher distances. As already men- 
tioned, however, the electron density can be probably retrieved, at best, only up to an unknown 
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Fig. 4. Same as Fig. 3, but using a model with both Bqi and C„ left as free parameters. 



multiplicative factor C„, representing an adjustment factor introduced to produce agreement with 
the actual RM observations (see also the discussion in Ingleby et al. 2007). Accordingly, in the sec- 
ond model (Case B), we performed a grid search by setting Bqi = (that is, assuming no dipolar 
component), allowing the parameter C„ to vary within a plausible interval of 0.3 < C„ < 3 at Steps 
of 0.01 Rq. This model was implemented to test the possibility that the radial component of the 
magnetic field might be already oc r^^ at distances as low as a few solar radii from the Sun (e.g., 
Vrsnak et al. 2004; Spangler 2005; Ingleby et al. 2007 ). The best-fit parameter to the observed data 
for Case B is C„ = 2.00 for the radial computation with R^s - 2.5 Rq. Again, the radial computa- 
tion with Rss = 3.25 Rq yields higher values, while the classic computation with Rss = 2.5 Rq 
yields a very high;^'^ value, hinting to low C„ values. The scatter plot of predicted versus observed 
RMs, with predictions obtained by using the best-fitting parameters for this case, assuming a PFSS 
radial computation with Rss = 2.5 Rq is shown in Fig. 3. Considerations similar to those concerning 
Fig. 2 are valid (see the discussion at the end of the previous subsection). 



4.3. Model with Bqi ^OandC„i^l (Case C) 

In the last model (Case C), we further performed a grid search allowing the two unknown parame- 
ters B()i and C„, to vary within the interval of < Bqi < 40, and 0.3 < C„ < 3. Also in this case, 
a minimum was very well defined allowing us to accurately determine the best-fit parameters. 
The model parameters that provide the best fit to the observed data are Bqi = 4.84 and C„ = 1.39 
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Table 1. The best-fit models. See text for more detail. 
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for = l/<T^^j-, and Bqi = 6.46 and C„ = 1.28 for a», = 1. Again, the best- fit was obtained for 
the radial computation with = 2.5 Rq. In Fig. 4, we show the scatter plot of predicted ver- 
sus observed RMs obtained by using the best-fitting parameters. Considerations similar to those 
concerning Fig. 2 apply even in this case. 

4.4. Discussion and comparison witti ottier models 

The model results are summarized in Table 1. For comparative purposes, for each considered 
model, the parameters C„ and Bqi that provided the best-fit to the observed RM's are tabulated 
together with the corresponding best- fit for both cases of non-uniform (a>, = l/f^jif , and a>, 
= 1 / cr* ^ .) and uniform (w, = 1 ) weighting, corresponding to the three steps outUned above for the 
evaluation of the uncertainties. We also report the estimate of a, for those models for which it is 
larger than 1%. We remark that the values of this variable strictly depend on the above hypotheses 
concerning the RM uncertainties (see the discussion before subsection 4.1), thus it has to be consid- 
ered with some care. For the model on the last line of Table 1 (classic computation with R^s = 2.5), 
we found that the is minimized by choosing very high values of the Bqi component, beyond 
the physical Umits expected to be plausible on the basis of present knowledge, so this model has 
probably to be ruled out. 

Overall, among the nine considered models (see Table 1), the three of them obtained by as- 
suming the position of the magnetic neutral line as deduced from the radial computation with 
Rss = 2.5 Rq always yielded the smallest minimum Xv for all weighting criteria. In case of the 
two weighting criteria of Steps 2 and 3, we obtained values ofxl at the minima small enough to 
allow for a meaningful statistical analysis (a > 1 - 2 %). The models computed by assuming the 
radial computation with Rss = 3.25 Rq also yielded well-defined minima within the allowed 
parameter space, although with larger x^ minimum values and with smaller a values with respect 
to the radial computation with R^s = 2.5 Rq. Finally, the classic computation always yielded the 
largest minimum (often at the border of the parameter space) and the worst or's among the three 
available WSO models. These results allow us to clearly draw our first and more robust conclusion. 
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that is, that the radial computation of the PFSS model from WSO with a source surface located at 
2.5 Rq is indeed the preferred choice, at least near solar minimum. 

Looking back at Table 1, it turns out that the minimum for Case C is always smaller than 
that for Case A and B for all weighting criteria. We remark, however, that this inequality is no 
longer true when comparing Case C and Case A by considering the xl^ instead of the due 
to the different v, although in both cases xl ~ ^ for the radial computation with R^s - 2.5 Rq. 
Related to these facts, the a value for Case A turns out, indeed, to be slightly higher than the 
one of Case C. However, the difference is quite small for both weighting criteria of Steps 2 and 
3, and thus we believe that both fits actually deserve full consideration. On the other hand, the a 
value of Case B appears slightly worse, especially in the case of the weighting criterium of Step 
3. Following the indications of our statistical analysis, we consider, for the sake of comparison of 
our results with other magnetic field estimates in hterature, both the models of Case A and C. On 
the other hand, taking into account the intrinsic uncertanties and assumptions related to the pB- 
inversion method used in this work, we are aware that our density estimate certainly needs some 
adjustment, at best by means of an enhancement or depletion factor. In fact, we believe that the 
slighdy larger a value for Case A with respect to Case C is just an artifact related to the fact that 
our x^ distributions are shallow. This impression is confirmed by comparison of our results in the 
two cases with independent work already published in the literature, as explained in the following. 

Our best-fit model for Case A, expressed as a power-law fit vahd in the range between about 5 
and 14 Rq, is 

B.= H|1.14^.7.40r-G, (5) 

r^ r'- 

whereas our best-fit model for Case C, expressed as a power-law fit vahd in the range between 
about 5 and 14 /?g, is 



4.84 1.43 



Br = ^ + ^^ 3.76 r""-"^ G. (6) 

Several empirical formulas have been introduced in literature in the past decades to model the 
radial magnetic field profile. In Fig. 5, we compare the magnetic field radial profile as derived 
from this work, given by eq. (5) and (6), with the radial profiles obtained in the inner heliosphere 
by Bulk & McLean (1978), Piitzold et al. (1987), Vrsnak et al. (2004), Gopaswamy & Yashiro 
(201 1) and Poomvises et al. (2012) with different techniques. Specifically, Dulk & McLean (1978) 
derived, from radio observations, a radial profile given by - 0.5 (r - 1)"'^ G for r < 10 Rq. 
The estimate by Patzold et al. (1987), B^ - 7.9 r"^^ G, valid in the range between 3 and 10 Rq, 
was also obtained through Faraday rotation measurements but using data from Helios. Vrsnak et 
al. (2004) used information on the band sphtting of type n radio bursts, obtaining Br = 1.4 r~^-^ 
G from the corona up to 1 AU. In the same plot, we also show further estimates of the magnetic 
field strength in the iimer heliosphere with various symbols (Sakurai & Spangler 1994; Spangler 
2005; Ingleby et al. 2007; Feng et al. 2011; You et al. 2012). From Fig. 5, it can be appreciated that 
the magnitude and variation trend of the radial component of the magnetic field as obtained in our 
analysis. Case C, is in fair agreement with the other estimates and in particular with the result from 
Patzold et al. (1987), at least in the range of heights around 5 Rq, being somewhat flatter overall. 
Moreover, the match with the estimate of Poomvises et al. (2012), obtained with a completely 
different method, i.e., from the analysis of the standoff distance of a CME-driven shock observed 
by the Solar Terrestrial Relations Observatory (STEREO) spacecraft, is particularly interesting. 
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Fig. 5. Comparison of the result from this work, eqs. (5) and (6), with various estimates (plotted 
with different symbols and colors) for the radial profile of the magnetic field strength in the inner 
heliosphere. 



On the other hand, the magnetic field profile we obtained instead by assuming C„ — I (Case A) 
appears to overestimate all other results, although its extrapolation at higher heights is still in fair 
agreement with the estimate of Poomvises et al. (2012). This substantiates our hypothesis that the 
He as estimated by means of the /jB-inversion method needs actually some correction (as suggested 
by the introduction of the enhancement factor C„), and that our best estimate of the radial magnetic 
field profile in the inner heliosphere is actually given by eq. (6) and not eq. (5). Finally, we remark 
that the radial profile obtained by Gopalswamy & Yashiro (201 1), who analyzed the shock stand-off 
distance and the radius of curvature of a flux rope associated to a CME event, is instead much flatter 
with respect to our profile. When comparing the above radial magnetic field estimates, however, we 
remind that all of them were obtained through remote sensing techniques and are therefore strongly 
influenced by the assumed (or, at best, estimated) electron density distribution along the LOS. 



5. Summary and conclusions 

Faraday rotation measures (RMs) estimated by Mancuso & Spangler (2000) along thirteen LOS to 
ten extragalactic radio sources occulted by the corona were used to constrain the inner heliospheric 
magnetic field around solar minimum. Since RM observations basically probe the electron density- 
weighted magnetic field strength along the integration path, the crucial point in this kind of analysis 
is obtaining a reliable estimate of the intervening electron density distribution along each LOS. 

By inverting LASCO/SOHO pB data taken during the days of observations, we were able 
to disentangle the two plasma properties that contribute to the observed RMs, thus allowing the 
radial component of the inner heliospheric magnetic field, assumed to have a simple analytical 
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form, to be uniquely determined. By comparing observed and model RM values, using a best- 
fitting procedure, we found that this profile can be nicely approximated by a power-law of the form 
B, - 3.76 G, in a range of heights spanning from about 5 to 14 Rq. In order to obtain the 

above estimate, an enhancement factor of C„ ~ 1.3 was required for the electron density rig as 
inferred by means of the pB-inversion technique. The magnitude and variation trend of the radial 
component of the magnetic field as obtained in our analysis is in fair agreement with previous 
estimates. Finally, our analysis suggests that the radial computation of the PFSS model from WSO 
with a source surface located at Rss - 2.5 is the preferred choice near solar minimum. 

Future direct measurements of the magnetic field from the magnetometers on board NASA's 
Solar Probe Plus mission are expected to further constrain the magnetic field radial profile in the 
range of heights investigated in this work. 
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